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Abstract 

This document contains a thorough introduction to the subject of 
conventional Quantum Computing, and the unconventional method 
known as Adiabatic Quantum Computation. Testing of a specifically 
written general simulation of the Adiabatic Quantum Computer with 
1, 2 and 3 qubits gives realistic results. Finally the application of an 
NP hard problem (Max Independent Set) to 4 and 8 qubit cases yields 
not only the cardinality of Max Independent Sets but also a visual 
representation. Leading to the development of a general algorithm for 
the Max Independent Set Problem that is executable on an Adiabatic 
Quantum Computer. 
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1 Introduction 



Adiabatic Quantum Computing (AQC) is a relatively new subject in the 
world of quantum computing, let alone Physics. Inspiration for this project 
has come from recent controversy around D-Wave Systems in British Columbia, 
Canada, who claim to have built a working AQC which is now commercially 
available and hope to be distributing a 1024 qubit chip by the end of 2008. 
Their 16 qubit chip was demonstrated online for the Supercomputing 2007 
conference within which a few small problems were solved; although the 
explanations that journalists and critics received were minimal and very lit- 
tle was divulged in the question and answer session. This 'unconvincing' 
demonstration^ has caused physicists and computer scientists to hit back at 
D-Wave. The aim of this project is to give an introduction to the historic 
advances in classical and quantum computing and to explore the methods of 
AQC. 

Initially Section [2] contains a short introduction to classical computation with 
very little concentration on the methods of implementing classical algorithms. 
The main topic in this section is the complexity theory used to classify prob- 
lems which in essence leads to problems that cannot be efficiently solved 
on classical computers, essentially justifying the need to develop quantum 
computers. 

Section [3] is a full introduction to the method of conventional quantum com- 
puting with explanations of basic principles including the DiVincenzo check- 
list and leading onto the implementation and failures of this method. Section 
[4] introduces the new unconventional method of quantum computing (AQC) 
with a thorough introduction of the problem that is extensively investigated 
in Section [U 

A brief summary of results with a final conclusion is available in Section 
[6} Suggestions of further research topics for future students are available in 
Section [7] and finally, the appendix containing only the program code of the 



1 http: / /www.scottaaronson.com/blog/?p=225 
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simulation used throughout this project is viewable in Section [8} 

2 Classical Computation 

Computers have come a long way since the very beginning. There have been 
countless physical realisations of the computer itself, moving from gears to 
relays to valves to transistors and finally to integrated circuits. All of these 
advances in technology have been extremely rapid and computers of this day 
and age are far more advanced than the early veterans of computing could 
have ever imagined. See Figure [T] for an eye-opening example. 

The perception of the slow and steady progression of computing evolution 
held until the publication of one paper in 1965, written by Intel co-founder 
Gordon E. Moore pQ. The general outcome of this paper was what is known 
now as Moore's Law which describes the trend in the history of computer 
hardware, that: 

'the number of transistors that can be inexpensively placed on an 
integrated circuit is doubling approximately every two years.' 

Although now the figure of 18 months is more widely accepted, the law looks 
set to continue for at least another decade. 

Even with the advances in technology and the continuation of Moore's Law, 
the only real change to computing is the greatly increased speed of compu- 
tation. The capability of the machines has yet to be altered. 

2.1 Complexity Classes 

Even with the greatly increased speed of computation, many problems are 
still unsolvable on the computers used today. For a computer to solve a 
problem, it must follow a specific set of instructions that can be computa- 
tionally applied to yield the solution to any given instance of that problem. 
This is known as an algorithm. The more efficient the algorithm the faster 
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'Scientists from the RAND corporation have created this model 
to illustrate how a "home computer" could look in the year 2004- 
However the needed technology will not be economically feasible 
for the average home. Also the scientists readily admit that the 
computer will require not yet invented technology to actually work, 
but 50 years from now scientific progress is expected to solve these 
problems. With teletype interface and the Fortran language, it will 
be easy to use.' 



Figure 1: Mock- Up of How a Computer Would Look in 2004 from 1954 
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the computer can solve that problem. Problems are categorised into classes 
of complexity, the class of complexity is determined by taking the most effi- 
cient known algorithm for a specific problem and calculating the number of 
computational steps needed to solve the problem as a function of the input 
size. Depending on how the number of computational steps relates to the 
size of the problem, the complexity class is found. This method is used, as to 
a computer scientist the crucial information about a problem is how quickly 
the time needed to solve it grows as the problem size increases. 

Below is a short introduction to various classical complexity classes, with 
official classifications taken from the National Institute of Standards and 
Technology 

NP Problems 

"The complexity class of decision problems for which answers can 
be checked by an algorithm whose run time is polynomial in the 
size of the input. Note that this does not require or imply that an 
answer can be found quickly, only that any claimed solution can 
be verified quickly. "NP" is the class that a Non-deterministic 
Turing machine accepts in Polynomial time." 

The class NP encapsulates all P and NP-Complete problems. The simple 
explanation for the problems contained in the NP class is that they can be 
checked in polynomial time. Meaning that all problems in the sub-classes of 
NP can also be checked in polynomial time. 

P Problems 

"The complexity class of languages that can be accepted by a 
deterministic Turing machine in polynomial time." 

Problems in the class 'P' are the type that a classical computer can solve in 
a 'reasonable' amount of time. The 'P' stands for polynomial time, which 

2 http://www. nist.gov 
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means that the time taken to solve the problem with the most efficient al- 
gorithm is the size of the problem raised to some fixed power. One of the 
simplest polynomial time algorithms is multiplication, where if we take two 
n-digit numbers the amount of time required to multiply those numbers is 
the total number of digits squared in 2 ). 

NP-Complete Problems 

"The complexity class of decision problems for which answers can 
be checked for correctness, given a certificate, by an algorithm 
whose run time is polynomial in the size of the input (that is, it 
is NP) and no other NP problem is more than a polynomial factor 
harder. Informally, a problem is NP-Complete if answers can be 
verified quickly, and a quick algorithm to solve this problem can 
be used to solve all other NP problems quickly." 

Problems in class NP-Complete are very different to NP problems, whereby 
every known algorithm for their solution will take an amount of time that 
increases exponentially with the problem size. Therefore solutions to these 
problems can be checked in polynomial time, but not necessarily solved in 
polynomial time. That is not to say that a polynomial time algorithm does 
not exist for these problems, it may yet be found, although for now this is 
the most widely accepted conjecture. The theory of NP completeness was 
developed by Cook, Karp and Levin in the 1970s [2], who have theorised that 
if an efficient algorithm for any one of these problems were to be found, it 
could be adapted to solve all other NP problems^- The implication of this 
is phenomenal, as developing a polynomial time algorithm for just one of 
these problems is enough to relieve many mathematicians of many years of 
anguish. 

A subset of NP-Complete is NP-Hard, these are NP-Complete problems in 
optimisation form. For example the travelling salesman problem is to find 
a route that visits all possible towns on a travelling salesman's list that in- 
duces a mileage less than some nominal value, which is NP-Complete. The 
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optimisation form of this problem is to find the shortest route possible whilst 
visiting all towns. 

One must stress that the explanations of complexity classes above are just 
conjecture. Results to date indicate that NP ^ P although this has never 
been proven true or false, therefore the question still remains open. In the fu- 
ture a classical polynomial time algorithm may be found for an NP-Complete 
problem, this would indicate that all NP problems are solvable in polynomial 
time. 

3 Quantum Computation 

The idea of quantum computing itself was loosely created by renowned physi- 
cist Richard Feynman in 1982(1]. His idea of simulating quantum mechanical 
problems with other more easily accessible quantum devices is the basis of this 
new branch of computing. Although initially these ideas did not spur much 
excitement between physicists in the form of actual quantum computation, 
although some research groups started investigating simple physical quantum 
systems that could possibly be used for this type of application. 

It wasn't until 1985 when David Deutsch published his theoretical paper on 
the universal quantum computer and what it might achieve that physicists 
started to take the subject more seriously. Deutsch indicated the possibility 
of simultaneous operations due to quantum superposition and the computa- 
tional efficiency of such a facility. He went on to develop his own algorithm, 
known as the Deutsch Algorithm that in simple terms can determine whether 
a (quantum) coin is bias with a single toss. Though it had little practical 
use, the race was now on to find a practical algorithm for use on a quantum 
computer. 

Again, interest in the subject soon dwindled until 1994 when Peter Shor de- 
vised an algorithm for efficient factorisation on a quantum computer that 
was published in 1996 [5]. This regenerated interest in the subject and ul- 
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timately led to the the next major discovery in 1997 by Lov Grover who 
devised an algorithm that can search a database containing N entries with 
approximately y/~N steps [BJ. 

Now thousands of computer scientists and physicists are attempting to de- 
velop new algorithms for use on physical quantum computers that are being 
developed by many companies all over the world. IBM are working on their 
NMR quantum computer, the Innsbruck group with their Ion-trap quantum 
computer, D-Wave systems with their superconducting qubits realisation, 
and many more. 

3.1 How Quantum Computing Differs 

New technologies will always allow classical computers to advance, although 
whether it is done by cramming more transistors onto a silicon chip or by 
increasing the clock-speed of processors the same fundamental fact that only 
the speed of computation is increasing still holds. The major fundamental 
difference with quantum computers is that the classical bit is replaced by the 
qubit (quantumbit) . This has created a variety of new possibilities. 

As a crude example, 4 classical bits can be re-arranged into 16 (2 n , n being 
the number of bits) different combinations (e.g. 0001, 0010, 0011, etc.). For a 
classical computer, these 4 bits combined can only be in one of the 16 possible 
states at any given time. The difference for a quantum computer containing 
4 qubits is that the qubits can be in a linear super-position off all sixteen 
possible states at any one time. Now say that both the classical and quantum 
computer perform a computation whereby all of these combinations must be 
explored. The classical computer will initialise the bits and perform the 
required operation to each individual state. Whereas the quantum computer 
can prepare the qubits in a linear superposition of all 16 states and perform 
the operation only once. 

At first glance this does not appear to be terribly impressive, until you re- 
alise the implications of how many more parallel operations can be achieved 
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in one clock cycle by simply adding more qubits. Just by adding a single 
additional qubit to any system, the number of possible states for that sys- 
tem doubles as does the number of possible simultaneous operations. For a 
classical computer to match this computational power, for every bit that is 
added, you would have to double the number of physical processors within 
the computer. 

3.2 Quantum Complexity Class 

The quantum complexity class is known as BQP (bounded error, quantum, 
polynomial time), this includes all problems that can be solved with a prob- 
ability of more than half in polynomial time on a quantum computer. Where 
BQP fits in with other complexity classes is best viewed pictorially in Fig- 
ure [2} It is widely believed that quantum computers will only prove useful 
if there are important problems in BQP that are not in P[7|. An example 
of a problem that has theoretically already entered this class is prime fac- 
torisation by Shor's Algorithm. Interestingly, this is probably the problem 
that is encouraging the development of quantum computers, as having the 
ability to perform fast prime factorisation will enable the deciphering of the 
RSA encryption used for the majority of monetary transactions across the 
internet. 



The qubit can be of many physical forms. Generally the qubit is a two state 
(or two states in a larger spectrum) quantum system, some examples include 
the vertical and horizontal polarisations of single photons, or the two spin 
states of a spin half particle. The state of which is described by \ip): 



Where 7, 9 (o<e<7r) and <fi (o<0< 2?r) are real. Although the overall phase factor 
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^ Computable 
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P solution can be found in time polynomial in the size (number of bits) 

NP solution can be checked in polynomial time 

NP COMPLETE any NP problem can be reduced to one of these 

BPP solution in polynomial time at probability p>1/2 

BQP solution in polynomial time at probability p>M2 on quantum computer 

PSPACE solution requires a polynomial amount of memory 

EXP solution can be found in exponential time 

COMPUTABLE solution can be found eventually 



Figure 2: Venn Diagram of Complexity Classes in Currently Believed Struc- 
ture 
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(e n ) is excluded as it has no observable effect [8]. Therefore Equation [T] can 
be re-written in its normalised form as follows: 

|V)=A|O)+0|1) (2) 

Where: 

|A| 2 + M 2 = 1 (3) 

This allows for a superposition of the two states |0) and |1), with |A| 2 and 
|/i| 2 as probabilities of measuring the classical states of and 1 respectively. 
The normalisation condition in Equation [3] permits the use of a geometric 
representation known as the Bloch sphere (see Figure [3]). This is a unit 
sphere where any point on the surface of the sphere is a possible pure state 
of the qubit it represents. 

As a qubit can only be measured in one of the two possible states, they are 
taken as givens. This has lead many to adopt the simpler vector notation for 
the state. The general state vector is given as: 

Where the top value is the linear coefficient for state |0) and the bottom |1). 
This representation will be used from now on. 

Another representation to consider is the density matrix. This is a far more 
general representation of the state of the system. The density matrix (p) can 
be used to represent mixed states (statistical mixtures of various states or a 
quantum system that has decohered) as appose to the ability of state vectors 
for sole representation of pure states only. The density matrix is formulated 
as follows: 



i 



(5) 
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Figure 3: The Bloch Sphere 
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Where is the state of qubit i. This representation will not be used for the 
majority of this document as the quantum computer system will be modelled 
without noise and decoherence for which the state vector representation will 
suffice. 



3.4 The DiVincenzo Criterion 



Many various two state systems could be classed as qubits, although for them 
to be useful in quantum computation other factors must be explored. The 
DiVincenzo criterion [9 J is a basic checklist outlining the properties that a 
quantum device must incorporate for it to be classed as a working quantum 
computer. The most basic checklist contains five criteria that are needed for 
quantum computing to be implemented although another two necessities are 
listed for the creation of an all-in-one quantum computing/communication 
device. 

1. A Scalable Physical System with Well Characterised Qubits 



The basic principles of the qubit have already been mentioned in Sec- 



tion [373] Generally if the qubit used satisfies the outlined description it 
can be classed as well characterised and may be useful for some imple- 
mentation of quantum computing. Although the scalability is one of 
the greatest issues that is yet to be addressed for many physical realisa- 
tions of a quantum computer, such as the Ion- Trap quantum computer, 
as currently this realisation has yet to reach double figures for the num- 
ber of coherent qubits. This requirement is needed to perform useful 
operations, as the number of accessible qubits limits the input size of 
the problem. 



2. The Ability to Initialise the State of the Qubit to a Simple State (e.g. 
|000...0» 
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In most cases the means of initialisation means cooling the system so 
that it is initialised in its known ground state, this must be achievable 
repeatedly with high fidelity, as from this state any other state can be 
reached by using the universal set of quantum gates (requirement 4). 
It is also possible to measure the qubits and once their state is known, 
operations can be performed to prepare them in their required initial- 
isation. This requirement is very important as without knowing the 
initial state of the system before computation, no sense can be made 
of an output after application of some algorithm. 

3. Adequately Long Decoherence Times 

The greatest worry for conventional quantum computers is decoher- 
ence. The largest source of decoherence is the entanglement of the 
system with the surrounding environment. These effects are detrimen- 
tal to the operation of a quantum computer. In many cases decoherence 
is described as a loss of information, although the information is only 
locally lost. In fact the quantum information stored in the system and 
environment is now shared, causing the qubits in a quantum computer 
to behave almost classically. In some cases decoherence has even been 
described as measurement of the system by the environment, although 
more worryingly it has even been shown that decoherence leads to an 
exponential increase of error rate with the input size[l]. 

It is widely believed that a sufficiently long time to perform a useful al- 
gorithm is > 10 4-5 clock cycles. Meaning that the quantum gates must 
be able to perform 10,000 operations within a single decoherence time. 
Many systems simply cannot achieve these expectations, although mea- 
sures have been made to increase the decoherence times of the quantum 
computer to an acceptable standard. It has been shown that quantum 
error correction codes can effectively be implemented [9]. Although sys- 
tems with error correction capabilities need a minimum of three qubits 
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to represent only a single one (to filter phase errors a minimum of 5 
qubits is necessary), this not only triples the number of qubits but 
significantly increases the number of resources needed to successfully 
complete an operation, including the physical 'gate' mechanisms. 

4. A "Universal" Set of Quantum Gates 

Quantum systems evolve unitarily according to their Hamiltonians, 
although for quantum computation to be performed on a specific sys- 
tem, the user must have full control of the qubits by manipulation of 
the applied Hamiltonians. This is due to the fact that the unitary op- 
erations build up the quantum algorithm, which are the specific set of 
instructions given to the qubits for any given problem before the result 
is read (just like classical computing). It is widely believed that any 
unitary transform can be formed from a series of single-spin operations 
and CNOT (controlled-NOT) gates [TDj. Therefore even if only these 
operations can be performed, this criterion will be satisfied. 

5. A Qubit- Specific Measurement Capability 

The final requirement is somewhat obvious. The need to be able to 
read ones answer (final state of the system) post computation is an 
ability that all quantum computers must possess. This requires a cou- 
pling between the quantum system and a classical measurement system, 
although this coupling cannot be continually apparent as this would in- 
troduce decoherence or even a complete collapse of the superposition 
state during the computation, deeming the computation a failure. 

Ideally the system will be in a pure superposition state after computa- 
tion, where all possible collapses of the state into a correct answer to 
the problem have a cumulative probability of (or very close to) 1. Yet 
this will not always be the case. It has been suggested that a trade-off 
between fidelity and other resources which result in reliable computa- 
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tion be made[9j. This is due to the simple fact that the target 100% 
efficiency will probably never be realised, therefore with systems that 
have a reliability/efficiency of 90% but need a higher figure (for exam- 
ple 97%) should simply be re-run to increase the accuracy of the results. 

To a lesser extent, the following two criteria would be useful for a universal 
quantum computing/communication device. If quantum computers ever be- 
come a household product, the dawn of a new era of hackers will be very real. 
With quantum computers the RSA encryption we all rely on in this day and 
age for monetary transactions across the internet will become easily decipher- 
able. For this reason, quantum key-distribution/cryptography will become 
necessary as will the need for a universal quantum computing/ communication 
device: 

1. The Ability To Interconvert Stationary and Flying Qubits 

2. The Ability to Transmit Flying Qubits Between Specified Locations 

3.5 Conventional Quantum Computing 

As with all computers, classical and quantum, the bits or qubits are initialised 
to some state. A series of logic gates are applied depending on the algorithm 
that the user wishes to implement and finally the output is read. Due to 
the quantum properties of the qubit, there are many more types of one or 
two qubit gates that can be implemented in a quantum system than possible 
in a classical system. These additional gates are designed to exploit the 
quantum properties of the system. For example, the Hadamard gate that 
is specifically designed to initialise a qubit into a superposition of its two 
classical states. 

Take a circuit diagram containing classical logic gates, to implement some 
algorithm on a desired input, where the input bit is represented by some 
voltage passing through a wire and the current state of this bit is dependent 
on the voltage. The bit of information travels through the wire and is acted 
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upon by the obstacles (logic gates) in its path. The circuit diagrams for a 
quantum algorithm acting on a selection of qubits look similar, although the 
qubits do not physically move throughout the system, the diagram is more 
of a timeline of operations that are applied to the qubits. Therefore a line 
on one of these diagrams will never make a U-turn for obvious reasons. 

There is currently no systematic way of programming a quantum computer 
and developing an algorithm for a specific task requires extreme 'cleverness', 
therefore there are very few useful algorithms that have been developed for 
implementation on a conventional quantum computer. The first and simplest 
quantum algorithm is known as the Deutsch Algorithm, it is used to check 
if a binary function is balanced with one evaluation of the function. If this 
were possible classically, we would be able to check if a coin was bias with a 
single toss. 

The algorithm is a series of gates applied to the qubits as pictured in Figure 
|1J The timeline is labelled with a series of unitary transformations (gates) 
that the qubits undergo, interestingly, all operations must be unitary so a 
quantum system is completely reversible unlike classical computers. The in- 
puts (initialised states) needed to implement this algorithm are two qubits 
each initialised as one of the two measurable states of a single qubit. There- 
fore l^o) — 1 1 ) , this state is then acted upon by two Hardamard gates which 
rotates both qubits into a superposition: 

l^i) is then fed into the two-bit controlled-f gate (Uf): 

, i 

|*2> = = - X>> ® (\f(x)) - \f(x))) 

x=0 

= ^E(- 1 ) /(x) i^)®(i°)-i 1 )) 

2 = 
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Figure 4: The Deutsch Algorithm in Circuit Diagram Form 



if /(0) = /(l) 
if /(0) 7^/(1) 




Now it is clear that \f(x))-\f(x)} = |0)-|1) if f{x) = and \f(x))-\f(x)) = 
|1) — |0) if f(x) = 1. Applying the Hadamard gate to the first bit again 
yields: 



|^ 3 ) = H®I\V 2 ) 




if /(0) = /(l) 
if /(0)^/(l) 



By measuring the first bit we can deduce that / is constant if the bit reads 
or balanced if the bit reads 1. 

As mentioned before, there are few useful algorithms that may be imple- 
mented on a physical conventional quantum computer. Most require far more 
qubits in the computation than the example above, and all involve many in- 
teresting sequences of gates to yield a final result. During a computation of 
this type, the applied gates transform the system into increasingly complex 
superposition states where a minimal disturbance can cause complete col- 
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lapse of the system and as mentioned in Section 3^ there is an exponential 
increase of error rate with input size. This has forced the development of 
new unconventional techniques. 



4 Adiabatic Quantum Computing (AQC) 

AQC differs greatly from the standard method of quantum computation. 
Rather than performing sequences of discrete unitary operations on one or 
two qubits individually over a period of time, a single time-dependant Hamil- 
tonian is applied to the whole system effectively reproducing the individual 
operations in a single evolution. Proofs of equivalence between the two meth- 
ods exist and more information on these can be found in Section 14.11 

To perform AQC, the system of qubits must evolve slowly in time from the 
ground state of an initial Hamiltonian with known eigenvalues, to a ground- 
state of the final Hamiltonian within which the answer to a specific problem 
is encoded, usually with unknown eigenvalues. There are many methods 
for this steady evolution although the most common is the slowly varying 
Hamiltonian [11] : 

H(t) = (1 - £) H I + TpHp (6) 

Where Hi is the initial Hamiltonian whose ground state is known (normally 
a linear superposition of all admissible states [12]) and Hp is the final Hamil- 
tonian within which the problem is encoded. Time t is the time during the 
evolution and T is the total time of the evolution. Larger values of T corre- 
spond to longer run times, and hence slower evolution. Thus the evolution 
should become increasingly adiabatic as T becomes large. It is clear that 
when t = the Hamiltonian is solely Hj and when t = T the Hamilto- 
nian is now Hp with smooth linear interpolation between Hj and Hp while 
< t < T. For simplicity the Hamiltonian is slightly altered, whereby t/T 
is replaced with fit where < t < 1 and the constant Q is a 'rate' that can 
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be changed for slower/faster evolution: 

H{t) = (1 - Sit) Hj + ntH F (7) 

The system of qubits evolves with time according to the Schrodinger Equation 
using the above time dependant Hamiltonian: 

-zh^(t)) = H(t)\^(t)} (8) 

By representing the state of the system with state vector notation and the 
Hamiltonian in matrix form, the solution to this equation can be found by 
solving the above as a system of ordinary differential equations (ODEs). This 
is the basis of the simulation described in Section 14.21 

AQC is thought to be a powerful method of solving NP-Hard minimisa- 
tion problems[l3]. The most common sector of mathematics to which these 
problems are related is Graph Theory, where nodes on a given graph are sub- 
jected to local biases as well as inter-node interactions. Therefore, in order for 
nodes to be represented by qubits there must be a method of implementing 
some local bias on the qubits as well as the inclusion of interactions between 
qubits. It is simply not feasible to experimentally implement problems of 
this type that involve fully connected graphs, therefore only nearest neigh- 
bour interactions can be modelled. The best known theory that describes 
these interactions is the Ising model[T3]. This model was initially developed 
to describe the inter-particle interactions in lattices to describe spontaneous 
magnetisation (amongst other effects). This model can be directly translated 
to be used on qubits as two-state magnetic spins of particles within lattices 
are directly analogous to the two states of any given qubit. 

Given a graph G, with vertices V and edges E the total Hamiltonian needed 
to describe this system is: 
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H F = J2 hoi 



(9) 



Where hi is the local Hamiltonian for qubit i and is the Ising interaction 
(coupling constant) between nearest neighbours i and j. Where if > it 
is energetically more favourable for qubits i and j to be aligned. Finally, a z is 
the Pauli spin matrix for the z-basis which is the classically measurable basis. 
More information on encoding specific problems into the form of the final 
Hamiltonian is available in Section 14.31 In most cases the local Hamiltonian 
is equal for all qubits, it may then be referred to as an external Hamiltonian 



In almost all cases the initial Hamiltonian Hj is taken to be of the form: 



As in all cases the ground state for this Hamiltonian is known or can be found 
easily. Where a x is the Pauli spin matrix in the x-basis. This initial Hamil- 
tonian is the reason why AQC is regarded as quantum, typically the ground 
state of qubits in this basis is a superposition state. This sets the method 
of AQC apart from classical annealing techniques which involve having an 
'always on' final Hamiltonian along with an external Hamiltonian of large 
magnitude in the same classical basis. Then slowly reducing the magnitude 
of the external Hamiltonian until only the problem Hamiltonian is applied. 
This method relies on the system switching between classical states to reduce 
the total energy of the system during the evolution until no further reduc- 
tion of energy can occur. After annealing, measurement of the system will 
frequently give the ground state, although in many cases the system remains 
in a local energy minimum state. Without raising the energy of the system 



(H): 




(10) 
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and neglecting possible quantum effects, there is no way of leaving a local 
energy minimum which is not the ground state of the system, and does not 
represent a solution to the encoded problem. Ruling this method of no use 
to any form of computation, simply due to its poor consistency. 

By starting in a ground state of the total Hamiltonian which is specifically 
constructed so that the ground state is a linear superposition of all possible 
states. Then slowly evolving the system so that it remains in a superposition 
of the possible ground states, local minima cannot exist. In the quantum 
method, states become less or more probable and do not switch between one 
another. One way of visualising this is to simply take the state vector of the 
system and observe the coefficients during a simulated evolution. They will 
simply increase or decrease depending on the Hamiltonian that is applied 
(although this is not experimentally possible). 

Assuming that a superposition state is easy to initialise and retain, the main 
source of error in an AQC is known as Landau-Zener tunnelling|15j. By 
evolving a system of qubits too quickly, it can tunnel into an energetic state 
that does not represent an answer to the encoded problem. The point in 
the evolution where this effect occurs most is at a so called 'anti-crossing', 
although this is in fact the point in the evolution where the eigenvalue gap 
between the ground state and the first excited state is at a minimum. A 
relationship to find an adequate evolution time from the minimum gap of 
the system has been directly derived from the adiabatic theorem [TT], and is 
as follows: 



Where £ is a constant for each system. Essentially, the size of T is governed 
by g^ n , although this is in no way helpful. With basic systems that can 
be effectively simulated by a classical computer, the eigenvalues can be pre- 
determined prior to simulation by checking all possible values for all possible 
times, although more complicated systems that are the sole reasons for the 
development of quantum computers will, in general, have unknown eigenval- 



T > 




nun 



(12) 
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ues prior to experimentation. Without prior knowledge of the eigenvalues, 
an adequate run time cannot be determined. 

Taking a single qubit evolution as an example, by evolving from a linear 
superposition of the two states (ground state of a Hamiltonian acting in 
the x-basis), to one of the two classical states (ground state of a Hamiltonian 
acting in either the positive or negative z-basis). Equation [7] for the evolution 
becomes a straight line equation through the Bloch sphere, the closer the 
state is evolved past the centre of the sphere (fully mixed state), the higher 
the possibility of tunnelling to an excited state. This means that errors will 
occur most if qubits are evolved between orthogonal states, as the eigenvalues 
cross during such evolutions. 

General decoherence that affects all quantum computers is due to entangle- 
ment with the environment. More information on the effect of this entan- 



glement can be found in Section 3.4 Although it seems that this effect is 



less detrimental to an AQC than a conventional quantum computer. Many 
studies and simulations of entanglement to the environment have shown that 
the AQC has an inherent fault tolerance and robustness against decoher- 
ence. The stray coupling to the environment affects the intended path of 
the eigenvalues of the Hamiltonian, although as long as a gap is present 
between the ground and first excited state the system will still evolve adi- 
abatically. In some cases this noise can even increase the fidelity of the 
system[T5] [16] [T7j [18] [19] [20] . For this reason, the effects of decoherence 
will be omitted from further study, and only perfect systems will be sim- 
ulated. 



With regards to the DiVincenzo checklist in Section 3^4, a system that can be 
used for AQC must as with conventional quantum computers, have scalable 
qubits, the ability to initialise and a measurement capability. As for the other 
criteria, the need for an adequately long decoherence time as explained above 
is not as necessary as it is for conventional quantum computers. Furthermore, 
the need for a universal set of quantum gates is specific to the operation 
of conventional quantum computers. Similarly, the ability to manufacture 
couplings between qubits and apply local Hamiltonians to individual qubits 
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is necessary for AQC. This ability allows the AQC to simulate any algorithm 
that can be applied to a conventional quantum computer which is equivalent 
to possessing a universal set of gates (see Section 4.1). 



4.1 Proof of Equivalence 

There is a somewhat simple proof of equivalence to conventional quantum 
computing by Mizel et al [21], although the interested reader may prefer the 
thorough proof by Aharonovj22] who has showed that any quantum circuit 
can be simulated by an adiabatic quantum algorithm. 

Interestingly these proofs show that an AQC is capable of performing all op- 
erations that are executable on a conventional quantum computer. Although 
no proof exists to show that a conventional quantum computer can repro- 
duce the same operations as an AQC. Due to this phenomenal turn of events, 
much interest has been generated among physicists and computer scientists 
for research into the AQC. To the extent that it is now believed that the AQC 
is a form of the universal quantum computer theorised by Deutsch[T6]. 



4.2 Simulating the AQC 

The simulation of the AQC has been written with the use of Matlab and the 
full code can be found in the Appendix. The simulation in basic terms solves 
the Schrodinger Equation as a system of ordinary differential equations. This 
is achieved by the user listing all Ising interactions between qubits and the 
local Hamiltonians on all the qubits as these are system dependent. The 
program will then compile the initial, problem and total Hamiltonians before 
finding and applying the initial ground state of the system. The program 
continues by compiling the ODEs and saving them to a separate '.m' file, 
these are then called and solved. The solver used is ODE45 which is based 
on an explicit Runge-Kutta formula known as the Dormand-Prince Pair, it is 
advised by Mathworks (creators of Matlab) that the ODE45 should always 
be used first and foremost. 
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The graphs that are plotted are dependent on the size of the inputted sys- 
tem. For all sized systems, the probability for individual states is plotted as 
a function of time, the eigenvalues of the system are plotted over time as well 
as an individual plot for the minimum gap energy between the two lowest 
eigenvalues. Furthermore, if the system is a single qubit, a plot of the evolu- 
tion of the state trajectory through the Bloch sphere is provided, and for a 
two qubit system the concurrence (measure of entanglement) is represented 
graphically as a function of time. Once simulation is complete there is an 
output of all states and their final probabilities in the Matlab window. 

4.3 Max Independent Set Problem 

A general algorithm for compiling a final Hamiltonian for any NP-Hard prob- 
lem to use on an AQC has yet to be developed or may not even exist. For 
this reason a specific problem that is known to be NP-Hard has been selected 
for further study. The chosen problem is known as the Max Independent Set 
(MIS) problem. Further constraints and adaptations have been introduced 
so that it is ideal for the simulation in this situation. At time of printing the 
independent set problem remains an NP-Complete problem for cubic pla- 
nar graphs[23], and its optimisation counterpart the MIS problem remains 
NP-Hard. 

The independent set problem is best described with graph theory. Taking a 
graph G=(V,E) with vertices (nodes) V and edges E, does an independent set 
exist that has a cardinality of at least k (k = some integer)? An independent 
set is a subset of the system of which no vertices are directly connected by 
an edge E. The most basic algorithm to solve this problem is to take every 
subset in the system that is of size at least k and examine each induvidually 
to determine if it is an independent set. The optimisation equivalent to this 
problem as mentioned above is the MIS problem, when the maximum k for 
the given graph G must be determined. This is the problem that will be 
applied to the AQC simulation. 
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Using the Ising model it is NP-Hard to calculate the ground state of a planar]^] 
two-dimensional system, whose inter qubit couplings favour anti-alignment 
of states whilst in the presence of an external Hamiltonian. If and only if 
each qubit is coupled to no more than three others [21]. At the same time it 
is NP-Hard to find the MIS of a cubicQ planar graph. Due to the equivalence 
of these problems the Ising model can be isomorphically mapped onto the 
MIS problem |25j. Whereby finding the ground state of the Ising model will 
also find an example of a MIS. For most problems there is a selection of 
MISs, but obtaining any one is NP-Hard[25j. On paper this ability yields 
nothing of interest, as both problems remain NP-Hard to solve. Although 
using a physical quantum system that behaves according to the Ising model, 
an instance of the MIS can be encoded into an applied Hamiltonian. This is 
useful as previously explained it is experimentally not feasible to implement 
a highly connected graph to any physical realisation of a quantum computer 
and secondly a planar graph can be implemented on a 2-dimensional chip 
containing qubits with only nearest neighbour interactions. Therefore, by 
measuring the ground state of the system, the MIS of thus system will also 
be determined. This is achieved by representing nodes on a graph with 
qubits, whereby the parity of the node is represented by qubit states |0) and 

|1>- 

Without mathematical proof, the mapping of the Ising model with couplings 
that favour ant i- alignment onto the MIS problem appears completely analo- 
gous (even without the presence of an external Hamiltonian) . If it is energet- 
ically favourable for qubits to be anti-aligned, this in almost all cases would 
cause every qubit to be of opposite parity to its nearest neighbours. If this 
were the case, finding the ground state of the system with the exclusion of 
an external Hamiltonian would still yield the MIS. Although this is not the 
case, many ground states of this form contain coupled pairs of qubits with 
the same parity that are opposite to other pairs (see Section 5.4 and Figure 



14). This in turn reveals no information about the MIS as no set of qubits 



3 A graph is planar if the vertices do not intersect when drawn in a single two- 
dimensional plane. 

4 The graph is known as cubic if all its vertices have three edges. 
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of the same parity are completely isolated from one another deeming them 
non-independent. For this reason, the local Hamiltonian must be introduced 
on all qubits, to act as some form of penalty function, so that only qubits of 
one orientation can be aligned with their nearest neighbours. 

By exploring the energies of every possible combination of qubits that a 
single qubit can be connected to, a small range of Hamiltonian magnitudes 
can be determined for which the ground state of the system may represent 
a MIS of the system. For a cubic graph, whereby each qubit may only have 
three nearest neighbours there are four possible combinations of states that 
the nearest neighbours may occupy (see Figures 5a to 5d where filled circles 
represent Is). By flipping the centre qubit between both states (0 and 1), 
and determining the total system energy of both orientations with Equation 
10] will allow inequalities to be developed. 



Assuming the qubit that is part of the independent set will be of opposite 
direction to the applied Hamiltonian, the state will represent independent 
qubits. The inequalities are determined by deciding which orientation of 
the qubit best describes an independent set, as the energy involving this 
orientation must be lowest of the two. Table [T] has a summary of the devel- 
oped inequalities and the Hamiltonian value that satisfy the inequalities if 
Jij = — 1 for all interactions. 

In summary 1 < hi < 3 for all i if = — 1 for all interactions. This also 
shows that no matter the size of the system, a Hamiltonian within this range 
will cause the ground state to be a MIS as it satisfies all possible interactions 
in that system. 



Situation 


Inequality 


Hamiltonian 


Fig. 


5a 




-AH - 3 J > -2H + 3 J 


H < 3 


Fig. 


5b 




J > -2H - J 


H > 1 


Fig. 


5c" 




2H - J > J 


H > -1 


Fig. 


5d 




AH — 3 J > 2H + 3J 


H > -3 



Table 1: Table of Inequalities to Determine a Hamiltonian Range 
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5 Results & Discussion 

Below is an extensive set of results for 1, 2, 3, 4 and 8 qubit cases. For one 
and two qubits, simple simulations to test the theory of adiabatic evolution 
as well as quantum properties of the system are explored. The three qubit 
case is used to investigate the effects of frustration in a quantum system. For 
4 and 8 qubits the MIS problem is applied, with exploration into the input 
parameters and the final results they yield. 

5.1 1 Qubit Simulation 

The system in this case has two possible classical states: 1 or 0. The initial 
Hamiltonian will be a linear superposition of the two states, this will be of 
the form: = + |0)). As there cannot be any Ising interactions, 

only a local Hamiltonian can be applied to the system. The Hamiltonian can 
work as a bias in either direction so these possibilities are explored as well as 
having no Hamiltonian whatsoever. The simulation is run for various values 
of Q of which the reciprocal is the number of units of time. 

It is clear from Table [2] that the simulation is providing the desired effects. 
When the local (problem) Hamiltonian is increased over time, the qubit aligns 
with the direction of the Hamiltonian as to remain in the ground state. This 
is also a good test of the adiabatic theorem; it is clear that with sufficiently 
slow evolutions the qubit is more likely to be measured in the ground state 
of the Hamiltonian at the end of the evolution. With a fast evolution that 
takes 0.1 units of time the qubit has an almost equal probability of being 
measured in either of the two states. Whereas after 10 units of time the 
likelihood of measuring the qubit in the ground state is greater with almost 
a probability of 1. This probability is further increased by running the sim- 
ulation for 20 units of time, although there is an interesting decrease after 
100 units. The combined probability of both states is therefore calculated 
for all simulations (with the values given in Table [2]). There seems to be a 
very small discrepancy of the total probability not accumulating to 1 after 
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Table 2: Table of Obtained Final Probabilities for Single Qubit States |0) 
and |1) in Specified Circumstances 



lengthier simulations. This effect appears to be dependent on time and the 
most legitimate reasoning is the rounding error introduced by the simula- 
tion's ODE solver. If a finite error is introduced into the probability after 
its calculation for every time interval, this error will become more prevalent 
with more time steps. Although in general the discrepancy is so small that no 
further action must be taken to obtain results of greater accuracy in future 
simulations. 

To obtain further understanding of the qubits reaction to the applied Hamil- 
tonians, plots of various aspects of the qubits evolution are produced for 
the application of the Hamiltonian that introduces a bias in the positive z- 
direction. This has not been done for the bias of -1 (as the plots are identical 
but with the states swapped). Nor have the plots been produced for ap- 
plication of no bias as this effect is only investigated to gather information 
on leakage and barely caused the qubit to evolve at all, graphically yielding 
nothing of interest. Plots are produced for the probabilities of states for 
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various values of Q (Figures 7a and 7b), eigenvalues of the two qubit states 



(Figure 6a) and the trajectory of the state of the qubit through the Bloch 



sphere as a function of time (Figures 8a and 8b). 



5.1.1 Larmor Precession 

By looking at the trajectories through the Bloch sphere for the single qubit 



case in Figures 8a and 8b and reviewing the plots of probability over time 
in Figures [7a] and 7b It is clear that some sort of oscillations are prominent 
in the probabilities, similarly there is a precession of the qubits orientation 
around the Hamiltonian direction in the Bloch spheres. This motion is similar 
to the trajectory of a point on the rim of a cone rolling on a flat surface. 
Further research reveals that the motion is due to a physical phenomena 
known as Larmor Precessions! 261. 



Take the qubit to be a spin-| magnetic moment where the Hamiltonian would 
be an applied magnetic field. If the magnetic moment is initialised at some 
direction and the field applied at an angle to this direction the field will 
exert some torque on the magnetic moment. Causing the moment to precess 
around the field direction, somewhat analogous to an ordinary compass with 
no oil (un-damped) in the earth's magnetic field or a spinning-top in the 
earth's gravitational field. 

To test the simulation further, a simple scenario can be run to directly view 
the Larmor Precessions. By taking a constant Hamiltonian in the z-direction 
and initialising the qubit to a state other than its ground state, a precession 
around the field direction along the edge of the Bloch sphere should be ob- 
served, mapping out a plane across the sphere perpendicular to the z-basis. 
The qubit must be initialised to a state other than that of the ground state 
as otherwise the precession would not be visible. The precession would be 
occurring although it would be aligned with the field, effectively turning on 
a spot. 

This raises the question as to why the precession is apparent during the initial 
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Figure 6: Single Qubit Example of Eigenvalues and Their Gap as a Function 
of Time During the Evolution from Hj to Hp 
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Figure 7: Single Qubit Probabilities of States as Functions of Time After 
Slowly Applying a Bias in the Positive Z-Direction 
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Trajectory through Bloch Sphere 




(a) Q =0.1 



Trajectory through Bloch Sphere 
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evolutions in Figures [8a] and [8b| as the qubit is initialised in the ground state. 
This is due to the evolution from an initial to final Hamiltonian, for short 
periods of time the qubit will not be in the ground state as the Hamiltonian 



is constantly varying. This is why the precession is more visible on Figure 8a 



compared to Figure 8b as the Hamiltonian is varying at a greater rate. 



This scenario can be solved easily by hand, here we have the Schrodinger 
Equation with the constant Hamiltonian applied to the general state vector 
notation. Where h has been normalised out: 





cos^j 
sin(|)e^ 



(13) 




smijje'"'' / ai \ sin^e^ 



(14) 

By solving the system as two ODEs, a solution to the top half of the equation 



is easily found and is shown in Equation [16} This means that the coefficient 
for state |0) is constant so the probability of measurement for both states 
remains constant as it solely relies of the value of 9. 



^cos(f) (15) 



.'. cos(|) = constant (16) 

Going on to solve the second ODE, knowing that cos(|) is constant the 
following can be determined: 



sin(f)e^ = i— sin(f)e^ (17) 
z at 
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(19) 



From Equation 19 it is clear that the phase is dependent on the time t. 



This time dependance is not visible in Figures [9a] and [9b] as the precession 
overlaps many times to trace a circle on the edge of the Bloch sphere. 

The Larmor frequency for full revolutions of the Bloch sphere is a known 
constant for many physically realisable qubits. This constant may be a good 
means of calibration for the unit of time in the simulation dependant on the 
type of qubit that is being modelled. 



5.2 2 Qubit Simulation 

The 2 qubit simulation is run to test various factors. Namely how the con- 
currence is affected by the applied final Hamiltonian (mainly the Ising in- 
teraction), a simple check of consistency of the calculated eigenvalues and 
finally if the final states have probabilities that resemble the expected most 
probable results. 

From studying Table [2] for the single qubit simulation, it appears that the 
most adequate single time scale to run these simulations is with an Q value 
of 0.1. This is the point that the probabilities give a clear indication of the 
correct results for the specific simulation and with this rate there is very little 
truncation error. 

Table [3] is a brief outline of the results obtained for Ising interactions -1, 
and 1 with and without an (identical) local Hamiltonian acting on each of 
the qubits. The table contains both the final eigenvalues and probabilities of 
each of the states. By calculating the ground state energies of all 4 states in 
each individual set of conditions using Equation [9j the eigenvalues have been 
confirmed. As for probabilities, the states with lowest eigenvalues (ground 
states) have the highest probabilities. In the cases where two or more states 
have the same lowest eigenvalue their probabilities are equal highest apart 
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Figure 9: State Precession Around Bloch Sphere with Application of a Con- 
stant Hamiltonian in the Z-Direction and Initialising the Qubit to the States 
Stated Above 
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|00) 


2 


0.000007 


1 


0.000316 


3 


0.000159 


|01) 





0.002724 


1 


0.000126 


-1 


0.248075 


|10) 





0.002724 


1 


0.000126 


-1 


0.248075 


111) 


-2 


0.993582 


-3 


0.998174 


-1 


0.503156 



Table 3: Table of Obtained Final Eigenvalues and Probabilities for All 2 
Qubit States in Specified Circumstances (Sl=0.1) 

from the case where Hi = 1 and J^ = — 1 where three states all have an 
identical final eigenvalue. In this case the probability is shared but not 
equally. As expected the two states related by symmetry (|01) and 1 10) ) 
have equal probability but the sum of their probabilities is equivalent to the 
single probability of state |11). This is unexpected but the most plausible 
reason for this is the fact that state |11) had the lowest eigenvalue during the 
majority of the evolution, at the end of the evolution when the two additional 
states joined the lowest energy group the possibility of tunnelling to one of 
these became less probable therefore the possibly of remaining in this state 
is higher. 

The concurrence is a measure of entanglement between the qubits, it ranges 
between (not entangled) and 1 (fully entangled). The concurrence is cal- 
culated over time from the state vector of the system, if the state vector has 
the following coefficients: 

I * > = C|oo> |00) + C|oi) |01) + C| 10> 1 10) + C|n> 1 11) (20) 
The entanglement is calculated in the following way: 
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Plot of concurrence during evolution 




Time 

(a) fn=0 



Plot of concurrence during evolution 




Time 

(b) hi=l 



Figure 10: Plots of Concurrence for Fixed Value of J=-l and Varying Local 
Hamiltonian (hi) 



5 RESULTS & DISCUSSION 



39 



x 1 o 15 Plot of concurrence during evolution 
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0123456789 10 

Time 

(c) hi=0 



x 10 5 Plot of concurrence during evolution 




0123456789 10 

Time 



(d) hi=l 



Figure 10: Plots of Concurrence for Fixed Value of J=0 and Varying Local 
Hamiltonian (hi) 
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Plot of concurrence during evolution 




Time 

(e) hi=0 



Plot of concurrence during evolution 
0.18 | 1 1 1 1 1 1 




Time 

(f) hi=i 



Figure 10: Plots of Concurrence for Fixed Value of J=l and Varying Local 
Hamiltonian (hi) 
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C = 2(| C |oo) | • \c\n) | - |c| i) | • |c|io> |) (21) 



The plots of concurrence during the evolutions can be found in Figures |10a 
to 



lOf When there is no interaction between the qubits (J=0), it is clear 
that they are completely isolated and work independently. The plots of 
concurrence for this situation in Figures 10c| and lOd are simply due to the 



truncation errors in the program, with the magnitudes of concurrence that 
are present, they can be regarded as not entangled at all. 



Studying Figures |10a and 10b where the interaction between the qubits 



means that it is more favourable for the qubits to be anti-aligned indicated 
that when no local Hamiltonian is applied the qubits only interact with each 



other. This is indicated in Figure 10a as the concurrence reaches almost 1 



which is the point of maximal entanglement. When a local Hamiltonian is 
applied, it is equally favourable for the qubits to be anti-aligned as well as 
in state |11) where both are aligned in the field direction. The concurrence 
reaches a value of 0.5 which can be said to be only half entangled. There is 
good reasoning for this, it is thought that if two qubits are entangled that 
by measuring one, information is gained of the state of the other. In the 
case where the qubits are maximally entangled (no local Hamiltonians), by 
measuring the first qubit it can be determined that the second will be of 
opposite parity. In the case of half entanglement, after measuring the first 
qubit, if the readout is 0, then the second qubit due to entanglement will 
be anti-aligned. But if the readout is 1, than it cannot be determined if the 
second qubit is in state or 1. In basic terms this means that only half 
the measurements on qubit one will reveal the state of qubit two, therefore 
half-entangled. 

In the scenario where it is more favourable for the qubits to be aligned due to 
their Ising interaction (J=l), without a local Hamiltonian the concurrence 
reaches a value of 1. Whereby if the first qubit is measured, the second 
qubit will almost always be in the same state, this is observable in Figure 
|10e More complications arise when the local Hamiltonian is applied to both 



5 RESULTS & DISCUSSION 



42 



qubits, the entanglement appears to rise at the same rate as before but 
suddenly drops off a third of the way through the simulation and ends with 
a very low concurrence of about 0.04. It looks as if the application of the 
local Hamiltonians breaks the entanglement of the qubits. From Table [3] 
it is clear that the most probable state is |11) although this appears to be 
fully determined by the local Hamiltonians and not the interaction between 
qubits. Therefore, by measuring the first qubit as 1, the second qubit will 
most probably also be one, but this will mainly be due to the applied bias 
and not the interaction between the two qubits. 

5.3 3 Qubit Simulation 

The three qubit simulation is the lowest possible system where frustration 
occurs. This is one of the many outcomes of the Ising spin model that can be 
applied to Spin Glasses (highly frustrated materials). If atoms in a lattice are 
arranged to be highly connected with nearest neighbours in some triangular 
form, it is easy to find the ground state of the system if the interactions 
between the atoms are ferromagnetic, as it is energetically favourable for all 
the atoms to have their magnetic moments aligned in the same direction. 
This ease is not reflected in the antiferromagnetic state. 

Taking a classical three particle example, where each particle is always con- 
nected to the other two, in a triangular arrangement: If anti- alignment is 
favourable, both other particles should have opposite spin to the first, al- 
though as they are also connected directly, do they accept being in the same 
alignment or does one of them flip? If so which one, and once flipped do the 
other two experience the same issue? The easy answer is that the system is 
frustrated, and that any arrangement of the three particles where one is the 
opposite parity to the other two can be found at any one time. The process 
of flipping orientation is completely random and the system will never be 
satisfied with the arrangement, so the probability of finding any one of these 
arrangements should all be the same. 

Replacing these particles with qubits, the probability of measuring each of 
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the states will be identical to the classical picture. Although when using 
qubits, due to their quantum nature, before measurement they will be in a 
superposition of all ground states and measurement will cause their wave- 
function to collapse into one of the classical states. This is clear from Table 
[4] as all states have an equal probability of being measured bar the two states 
where all qubits have identical orientations. 

After applying the local Hamiltonian to all the qubits, three states become 
most favourable (|011), 1 101) and 1 1 10) ) . All of them include two Is and 
one and can be regarded as identical due to symmetry. If the triangle of 
qubits were to be rotated you would not be able to distinguish one state 
from another. The reason for these states having the lowest energies is due 
to the fact that having two Is lowers the energy of the system as two of the 
qubits are aligned with the local Hamiltonian direction which is energetically 
more favourable. Even though the third qubit is anti-aligned to the local 
Hamiltonian and raises the energy of the system, by being of this orientation 
more of the Ising interactions are satisfied, thus lowering the energy overall. 
The eigenvalues and probabilities of states listed in Table [4] can be found 
from the plots in Figures lla|to 12b 



5.4 4 Qubit Simulation 

The 4 qubit simulation is the lowest possible sized system for the application 
of a cubic planar graph. All qubits are linked to one another as shown in 



Figure 13 Meaning that this system is also the lowest sized system for 
the application of the specifically manufactured MIS problem. This is an 
interesting system to initially test as the ground state of this system without 
the application of an external Hamiltonian may in fact contain coupled pairs 



of qubits with the same parity that are opposite to other pairs (Figure 14). 
This in turn reveals no information about the MIS as no set of qubits of the 
same parity are completely isolated from one another deeming them non- 
independent. For this reason, the local Hamiltonian must be introduced on 
all qubits, to act as some form of penalty function, so that only qubits of one 
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(a) Eigenvalues 



Probabilities 

0.18 r 




0123456789 10 

Time 

(b) Probabilities 

Figure 11: Plots of Probabilities and Eigenvalues as Functions of Time for 
a Three Qubit System with Antiferromagnetic Interactions (f2=0.1, Jij=-l, 
hi=0) 
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(a) Eigenvalues 



Probabilities 

0.35 r 




Time 

(b) Probabilities 

Figure 12: Plots of Probabilities and Eigenvalues as Functions of Time for 
a Three Qubit System with Antiferromagnetic Interactions. (f2=0.1, Jij=-l, 
hi=l) 
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Oil) 
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-1 
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0.330406 
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0.000082 





0.007053 



Table 4: Table of Obtained Final Eigenvalues and Probabilities for All Three 
Qubit States with Antiferromagnetic Interactions (fi=0.1) 



orientation can be aligned with their nearest neighbours. 

Table [5] contains final probabilities of all admissible states in the absence 
of an external Hamiltonian acting on all qubits. The ground states of the 
system with no external Hamiltonian contains two Is and two 0s arranged in 
six possible combinations. These are clearly not independent sets as neither 
of the two 'sets' are independent. It is simple to deduce that the independent 
set for this 4 qubit system only contains one qubit that is anti-aligned to the 
other three, where the qubit that is anti-aligned is of opposite parity to the 
applied Hamiltonians direction to minimise the energy of the system. After 
introducing the external Hamiltonian (local Hamiltonian acting equally on 
all qubits) the results obtained are available in Table 16} 



As mentioned above, the independent set for this system has a total cardinal- 
ity of 1. Meaning that the independent set comprises only of a single qubit. 
It is known that the qubit will be of opposite parity to all other qubits as 
well as the local Hamiltonian that it is subjected to, indicating states that 
represent independent sets comprise of three Is and one 0. Where the qubit 
of state 0, represents the MIS. From the collected data in Table [6] it is clear 
that the states that represent the correct independent sets in fact have the 
highest final probabilities although share their eigenvalue with the states that 



have two qubits of each orientation (Figure 18a). This means that four cor- 



rect states and six incorrect states have the identical ground state energy. 
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Figure 13: Fully Connected Cubic Planar 4 Node Graph G=(V,E), where 
Vertices V Represent Qubits and Edges E Ising Interactions 




3 2 



Figure 14: An Example of a Ground State in a 4 Qubit System that Does 
Not Represent a Max Independent Set 
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Table 5: Table of 
Qubit States in the 



Obtained Final 
Absence of an '. 



. Eigenvalues and Probabilities for All 4 
External Hamiltonian 
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Table 6: Table of Obtained Final 
Qubit States in the Presence of an 



Eigenvalues and Probabilities for All 4 
External Hamiltonian 
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Using the results where £7=0.1, the probability of obtaining a correct answer 
is just above 0.5. To increase this probability, the parameters within the 
simulation must be adjusted. Measures must be taken to reduce the energy 
of the correct states and increase the energy of the incorrect states. Table [7] 
has a list of equations that calculate the energy eigenvalues of each type of 
state. 

In Table [7] the top two combinations of states have been separated with a 
horizontal line as no matter how much the magnitudes of the variables hi 
and Jij are increased, their total energy will only increase. Leaving only 
three other combinations of states. To lower the eigenvalue of the needed 
combination (three Is and one 0) to be less than the conflicting combination 
(two Is and two 0s) it is clear that the local Hamiltonian acting on the qubits 
must be increased. Although increasing it too much will cause the state 
1 11 11) to become more energetically favourable as the Hamiltonian strength 
will overwhelmingly supercede the interaction strength and cause all qubits 
to align in the Hamiltonians orientation. Traditionally the magnitude of J is 
equal to 1 (and negative for antiferrromagnetic interactions), continuing this 
trend only the parameter of h must be adjusted. Plotting the eigenvalues 
that may be obtained against values of h gives Figure [l5j Clearly the correct 
states are more energetically favourable when: 1 < hi < 3, as they have the 
lowest eigenvalue between this range. 

Knowing how the system will react to Hamiltonians of various magnitude, 
it is clear that further simulations are needed to investigate the probability 
of obtaining a state that represents a MIS for varying strengths of local 



State 


Eigenvalue 


|0000) (xl) 


4h- 6J 


0001) (x4) 


2h 


|0011) (x6) 


-2J 


0111) (x4) 


-2h 


1111) (xl) 


-4h - 6J 



Table 7: Table of Equations for the Calculation of Eigenvalues of Various 
States (J=-l) 
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Hamiltonian Magnitude 



Figure 15: Eigenvalues as Functions of the Magnitude of the Applied Local 
Hamiltonians (J=l) 



Hamiltonian and for a range of simulated times. The obtained data is plotted 



in Figures 16 and 17 



From Figures 16 and 17 it can be deduced that the adiabatic theorem is 
still a major factor in determining the state that the system evolves to. 
In short time periods the probability of finding a correct state is almost 
constant across all Hamiltonian magnitudes, with almost no dependance on 
the magnitude of the external Hamiltonian. This is also deducible from 
Table |8j as for a time of 0.1, the probability varies very little from 0.25 
for obtaining a correct final state under all tested Hamiltonian magnitudes. 
Although for this system out of sixteen possible states there are four that 
can be regarded as correct, exactly a quarter of the total number. Meaning 
that within this time period, all states right and wrong have almost identical 
final probabilities. 

As expected for Hamiltonian magnitudes between one and three the probabil- 
ities of measuring a state that represents a MIS becomes almost one for large 
simulation times (from 20 units of time and greater). With these magnitudes 
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0.3465 


10 


0.5184 


0.9396 


0.9816 


0.8836 


0.5704 


0.0944 


20 


0.474 


0.9648 


0.9924 


0.9912 


0.4604 


0.0348 


30 


0.512 


0.982 


0.9932 


0.988 


0.5592 





40 


0.482 


0.99 


0.9928 


0.986 


0.4716 


0.006 


50 


0.5016 


0.9916 


0.9916 


0.9884 


0.5308 


0.0048 



Table 8: Combined Final Probability of Correct Four States for Varying 
Time and Hamiltonian Magnitude 




Time 



Figure 16: Combined Probabilities of Independent Sets as Functions of Time 
for Varying Magnitudes of Local Hamiltonians (J=-l) 
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Figure 17: Probabilities for Obtaining an Independent Set for Varying Times 
and Hamiltonians in Contour Form (J=-l) 



the eigenvalues of the correct states split from the others with a reasonable 
gap. The states that start with the four lowest eigenvalues all converge to a 
single final eigenvalue meaning that the final state of the system is a linear 
superposition of these four states. The changes in eigenvalue progression are 
very interesting for the various Hamiltonian magnitudes, these are viewable 



in Figures 18a to 18f 



Figure [18a] is the situation with hi = 1 for all i. It is clear that a large set 
of states both correct and incorrect all evolve to become energetically equal. 
At the end of the simulation where all states converge, the probability of en- 
tering any of these states should be almost equal, although due to symmetry, 
the incorrect states have less probability individually. This is because the 
total probability of the incorrect independent sets (ground states) is equal 
to the total probability of the correct independent sets for sufficiently long 
simulation times (see Table [6l). 



Continuing to Figure |18b| where hi = 1.5 for all i the eigenvalues for the 
states that do not represent independent sets have split from those that do 
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Eigenvalues 




(a) hi=l 



Eigenvalues 




(b) hi=1.5 



Figure 18: Plots of Eigenvalues for Fixed Value of J=-l and Varying local 
Hamiltonian (hi) 
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Time 

(c) hi=2 



Eigenvalues 




Tims 

(d) hi=2.5 



Figure 18: Plots of Eigenvalues for Fixed Value of J=-l and Varying Local 
Hamiltonian (hi) 
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Figure 18: Plots of Eigenvalues for Fixed Value of J=-l and Varying Local 
Hamiltonian (hi) 
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by only one energy unit. Meaning that the probability of tunnelling into more 
excited states will be far smaller for sufficiently slow simulations. Similarly 



in Figure |18c[ where hi = 2 for all i the eigenvalues are still separated but 
with a larger gap, yet increasing the local Hamiltonian further to 2.5 (Figure 



18d) the gap reduces again to only one unit of energy. Meaning there is 



symmetry in the eigenvalue gaps in this system due to the magnitude of the 
Hamiltonians about a centre which is when hi = 2 for all i. Yet even with a 
smaller gap, it seems that there is very little difference between probabilities 



obtained using hi = 1.5, 2 and 2.5 (see Figure 16). Finally, Figure 18f is the 
plot of eigenvalues when hi = 3.5 for all i, in this situation there is only a 
single ground state, with a well-defined gap between the ground state and 
all other states throughout the evolution. This ground state as expected is 
the state where all qubits are aligned simultaneously with the Hamiltonian 
direction, as the lowering in energy due to the Hamiltonian magnitude is 
far greater in this situation than that of the Ising interactions. Therefore 
the probability of measuring a state that represents a MIS at the end of the 
simulation is dramatically reduced, as this relies on the system randomly 
tunnelling into a correct state. 



5.5 8 Qubit Simulation 
5.5.1 8 Qubit Cube 

The 8 qubit cube is studied as it is the simplest 8 qubit system due to various 
symmetries. The cube can be drawn in cubic planar form by simply taking 
a cube, which by default each qubit is linked to three others, and squashing 
it into a 2 dimensional plane. The qubits are situated at the corners of the 



cube as shown in Figure 19 



Every link in the system can be satisfied with antiferromagnetic interactions. 
Every qubit can be of opposite orientation to all of its nearest neighbours. 
Meaning that exactly half of the qubit will be of each orientation. An example 



of the exact MIS is viewable in Figure 20 where filled qubits are of one 
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Figure 19: 8 Qubit Cube in Cubic Planar 8 Node Graph G=(V,E) Form, 
where Vertices V Represent Qubits and Edges E Ising Interactions 




Figure 20: Example of the Max Independent Set of a Cube 
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orientation and the unfilled are the other. It is not important which ones are 
which due to the symmetries, by just rotating the system about its centre by 
ninety degrees the filled and unfilled orientations swap. 

As all of the ferromagnetic Ising interactions in this system can be satisfied 
simultaneously, no local Hamiltonian needs to be applied to the system in 
order for the ground state of this system to be representative of a MIS. 
Although this is a special case and not representative of all cubic planar 
systems, therefore this system is investigated to find a threshold Hamiltonian 
under which the system's ground state no longer represents a solution to the 
MIS problem. 

As mentioned before, due to symmetries there are only two degenerate ground 
states for this system, without an application of a local Hamiltonian on each 
qubit. These are |011001 10) and the opposite parity 1 10011001). To reach 
these, the system would be prepared in a superposition of all admissible states 
and the evolution would involve slowly increasing the Ising interactions whilst 
decreasing the magnitude of the initial Hamiltonian. The system by the end 
of the evolution would be a superposition of the two ground states (men- 
tioned above) and according to the Copenhagen interpretation of quantum 
mechanics, measurement would collapse the system's wavefunction into one 
of these states. As the two degenerate ground states are of completely op- 
posite parity, it would appear that the system is in a fully entangled state, 
so by measuring only one of the qubits the state of all others could be de- 
termined (if the ground state is known prior to evolution). Even though this 
entanglement is obvious from the two admissible ground states, the concur- 
rence of a system greater than 2 qubits cannot practically be determined |2 7] . 
Obtaining one of these states would obviously depend on the sufficiently slow 
adiabatic evolution and the robustness of the system to thermal excitation 
into a more energetic state. 

The obtained probabilities for various Hamiltonian magnitudes and simula- 
tion times are available in Table [9] By raising the Hamiltonian magnitude 
to three, the most probable state changes to the state where all qubits are 
aligned with their local Hamiltonian direction, which is the same for all 
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Hamiltonian Magnitude 


Time 





1 


2 


3 


0.1 


0.0082 


0.0082 


0.0082 


0.0082 


1 


0.0608 


0.0522 


0.0342 


0.0182 


2 


0.233 


0.2024 


0.0958 


0.0144 


10 


0.8374 


0.8496 


0.8704 


0.0186 


20 


0.9678 


0.9698 


0.9656 


0.0104 


33.33 


0.9822 


0.9822 


0.9822 


0.0322 



Table 9: Combined Final Probability of Correct Two States for Varying Time 
and Hamiltonian Magnitude for a Cube of Qubits 




Figure 21: Total Probability of Obtaining an State Representing the Max 
Independent Set After Measurement (J=-l) 
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qubits. Yet there is very little to distinguish between the probabilities ob- 
tained for identical simulation times for Hamiltonian magnitudes from zero 



to three, as seen in Figure 21 



Plots of the eigenvalues and probabilities of states for all of the four situations 



with varying local Hamiltonian magnitudes can be found in Figures 22a to 



22h It is obvious that the minimum gap between the two lowest eigenvalues 



will always be zero due to there being a pair of degenerate ground states, 



deeming Equation 12 for the relationship between minimum gap and evo- 
lution time of even less use. The equation suggests that the run time for a 
system with degenerate ground states should be infinite. There appears to be 
no criteria on using this equation when the degenerate ground states individ- 
ually represent correct answers to the encoded problem. If the gap between 
the ground state eigenvalue and first incorrect state eigenvalue is considered, 
there is more sense to the equation. Exciting the system to the first incorrect 
state will induce errors in the final measured value. The argument against 
such consideration is that one of the final degenerate ground states starts the 
evolution with its energy equal to a selection of incorrect states, where there 
is no energy gap. Although this would not affect the situation as at the start 
of the evolution the system is initialised to a non-degenerate ground state 
that is not energetically equivalent to an incorrect state. During the evolu- 
tion, if the system were to be excited into the state that eventually evolves 
to become a degenerate ground state, this would increase the probability of 
measuring a state that represents a correct answer. To enter an incorrect 
state, the system would have to be raised from the initialised ground state 
to the first excited state and again to the second excited state, which is en- 
ergetically equivalent to the difference between the initialised ground state 



and the second excited state. By considering Equation 12 with the gap be- 
tween initialised ground state and lowest incorrect state, it is clear why the 
probabilities for identical simulation times but varying local Hamiltonians 
are similar in Table |9j By observing the eigenvalue evolutions, the new sug- 
gested minimum gap is 2 (and not 0) for all local Hamiltonian magnitudes 
and this is not due to the magnitude of the local Hamiltonian applied to 
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all qubits, but a consequence of the initial Hamiltonian. Meaning that the 
minimum eigenvalue gap across all simulations in Table [9] is constant, giv- 
ing good justification for the almost identical probabilities obtained across 
various Hamiltonian magnitudes. 



5.5.2 Random 8 Qubit 



The random 8 qubit simulation in Figure 23 has been specifically devised 
and investigated to confirm previously established relationships in Sections 



5.4 and 5.5.1 By using the brute force method to find all possible MISs for 
this system, six different possibilities were found. Many related by symmetry 
as shown in Figures 24a and |24b 



As there are a possible 256 (2 n ) states that this system can occupy, so only 



the probabilities of the six MISs have been recorded in Table [TOj where Pi 
is the individual final probability of the two states in set 1 and P 2 is the 
individual final probabilities of the four states in set 2. P? is the total final 
probability of all 6 states that represent MISs. 



From Table [10] it is clear that the simulation has found the correct MISs, 
in some cases with a total probability of measurement of 0.9848 for 6 states 
out of a possible 256. This is good reason to believe that the simulation has 
worked, considering that if all states in this system were equally probable, 
the total probability of 6 states would amount to a value of 0.0234. 

Plots of the individual set probabilities and total probabilities have been 
compiled as functions of time and local Hamiltonian magnitudes. These are 



viewable as Figures 25a and 25b The plots of probability as functions of 



total time of the evolution in Figure 25a reveal that the adiabatic theorem 



still holds strong. Similarly to simulation in previous systems, the probability 
saturates to a value very close to 1 at large units of time, around 10 to 15 
units. Interestingly the individual sets appear to have more symmetry than 
initially expected. The two MISs in set 1 converge to have a probability 
of about 0.25 each at saturation and ideally in total this would amount to 
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Figure 22: Plots of Eigenvalues and Probabilities Whilst J=-l, hi=0 and 
Q = 0.1 (Legends Have Been Omitted as they Include 256 Entries) 
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Figure 22: Plots of Eigenvalues and Probabilities Whilst J=-l, hi=l and 
Q = 0.1 (Legends Have Been Omitted as they Include 256 Entries) 
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Figure 22: Plots of Eigenvalues and Probabilities Whilst J=-l, hi=2 and 
Q = 0.1 (Legends Have Been Omitted as they Include 256 Entries) 
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Figure 22: Plots of Eigenvalues and Probabilities Whilst J=-l, hi=3 and 
Q = 0.1 (Legends Have Been Omitted as they Include 256 Entries) 
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Table 10: Table of Obtained Final Probabilities for Independent Set States 
in Specified Circumstances 
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Figure 23: Random 8 Qubit Cubic Planar 8 Node Graph G=(V,E) Form, 
Where Vertices V Represent Qubits and Edges E Ising Interactions 

exactly half of the available probability of states. The four MISs in set 2 
appear to converge to a value of 0.125 each and their cumulative total would 
also amount to a value of 0.5 in an ideal system. By the end of the evolution 
all of the states that represent MISs have an identical eigenvalue of -10 and 
the states would be in some superposition with entanglement throughout the 
system. This entanglement is the most likely the cause for the interesting 
symmetry between states, although as previously explained, no method yet 
exists to calculate the entanglement between more than two qubits, and 
leaves this phenomena as a matter for further investigation. 



Figure 25b| shows the symmetry across various magnitudes of local Hamilto- 



nian. There are two effects occurring within this plot, firstly non-MISs be- 
coming more favourable and secondly the adiabatic theorem. When all lines 
converge either side of the centre near Hamiltonian values 0.75 and 3, these 
are the points at which other states (non-MISs) become more favourable. 
Interestingly at this point, the values for probability at T=l are the highest, 
this is induced due to the adiabatic theorem. With these local Hamiltoni- 
ans (0.75 > H > 3), other states become more energetically favourable, so 
evolving the system slowly the correct MISs measurement becomes less prob- 
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Figure 24: Graphs of All Possible Ground States Separated Into Two Sets 
By Symmetry 
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Hamiltonian Magnitude 



(b) Probabilities as Functions of Final Hamiltonian Magnitude 
Figure 25: Plots of the Probabilities Obtained for the Random 8 Qubit Graph 
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able, although evolving the system too quickly will cause almost all states to 
have an almost equal probability of measurement. Towards the centre of the 
plot, there is clearly some symmetry across the values of local Hamiltonian, 
although by simply evolving the system even slower will cause almost all 
probabilities to become close to 1, for a large range of Hamiltonian values, 
at a minimum between 1.5 and 2. 



6 Conclusion 

This project has exploited the isomorphic mapping of the ground states of 
multi-qubit systems (with anti-ferromagnetic coupling) to the well known 
MIS problem. This permitted an investigation into the magnitude of the 
required local Hamiltonian for the development of a general algorithm. It is 
clear that some energy penalty must be introduced via the local Hamiltonian 
to all qubits, so states that involve qubit pairs of the same orientation cannot 
appear in both parities. As in these cases neither parity can represent an 
independent set. 

Significant progress has been made on understanding the basic nature of the 
MIS problem. By observing the eigenvalue evolutions and the probabilities of 
measuring the ground state of the system under varying magnitudes of local 
Hamiltonian, a better assumption of the Hamiltonian that must be applied 
to yield a ground state that represents a MIS on all possible new systems 



can be made. For the 4 qubit system in Section |5.4[ it was found that a 
value for the local Hamiltonian between 1 and 3 was appropriate and values 
less than 1 or greater than 3 cause the ground state of the system to be a 
state other than one that represents a MIS. Furthermore, it was found that 
at Hamiltonian magnitudes of exactly 1 and 3, the various possible states 
that represent MISs shared their eigenvalue with states that did not, there- 
fore values between 1.5 and 2.5 are considered to be more acceptable. The 
simplest 8 qubit example that is referred to as a cube in Section |5.5.1| was a 



special case, where without an acting Hamiltonian the system has a perfect 
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ground state where all Ising interactions can be satisfied. The use of this 
system, introduced the possibility of applying a Hamiltonian where it is not 
necessary for the ground state to represent a MIS, this is useful to find a 
maximum threshold Hamiltonian. Essentially the point at which the energy 
penalty introduced by the local Hamiltonian (acting on all qubits) causes 
a ground state that does not represent a MIS to become most probable at 
measurement after the evolution. The threshold was found to be when the 
Hamiltonian magnitude was 3. This is good evidence that applying a Hamil- 
tonian to find a MIS can work with any system even if not entirely necessary, 
allowing a general algorithm to be established for this problem. The random 



8 qubit system in Section 5.5.2 simply confirmed previous results, although 



due to the sharpening of the peak by the adiabatic theorem in Figure 25b 



it appears that the symmetry in local Hamiltonian is not exactly centred 
about the value of 2. A basic Figure that can be read from the plot is about 
1.75. A solid relationship has yet to be established for the ratio of coupling 
strength to local Hamiltonian magnitude, although if the system is evolved 
sufficiently slowly a range of values can be used, which in turn can yield some 



very interesting results (see Section 7.3). 

The scalability of this problem to larger systems is expected to work from 



results obtained in Section |4.3[ although due to computational restrictions 
within the student version of Matlab no system over 8 qubits can be simulated 
with the general program that has been written for this project. Implemen- 
tation of the general algorithm on larger systems should not cause any issues, 
this is mainly due to the restrictions in place on the type of problem that can 
be simulated. For instance, the fact that the systems have been restricted to 
be cubic. Meaning that every qubit can only be linked to three others, this in 
turn allows for the scalability of the algorithm to larger systems, as only the 
three nearest neighbour interactions of every qubit need to be satisfied. 

Very little of this project has been spent on the time necessary to evolve a 
system adiabatically. The theorem has been shown to hold, as in all cases the 
results became increasingly accurate with larger evolution time. There are 
various reasons for not spending much time on this subject, mainly the fact 
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that the simulation remains to be calibrated, meaning that if any analysis 
were to be made, there could be no comparison made to any physically 
realisable system of qubits. Secondly, the main research topic of this paper 
was to investigate the implementation of an NP-Hard problem to AQC with 
the hope of developing a general algorithm for a specific problem. 

Finally, the MIS problem is the minimisation version of the NP-Complete 
problem known as the Independent Set problem. Although, this pair of 
counterpart problems is very unique for a single reason; solving the NP- 
Hard problem and finding the cardinality of the MIS for any system will also 
solve all instances of the NP-Complete counterpart. By simply determining 
the max k value of a graph G=(V,E), determining if k is greater of less than 
some nominal value is trivial. Hence solving all instances of the NP-Complete 
counterpart. 

In conclusion, to obtain a visual representation and the cardinality of a 
MIS of any given cubic planar system, one must apply a final Hamiltonian 
in an adiabatic evolution that contains only interactions that favour anti- 
alignement (Jjj = —1) with an applied Hamiltonian of magnitude 2 (hi = 2) 
for all qubits. This is in agreement with the mathematical analysis in Section 

7 Recommendations for Further Research 

There are many areas of AQC that remain to be explored. Below is an 
extensive list of various effects that could in theory either quicken or improve 
the results of AQC. 

7.1 Quantum Annealing 



With reference to the 8 qubit cube simulation in Section 5.5.1 the proba- 
bilities gained by the standard method are expected to increase by doubling 
the magnitude of the initial Hamiltonian, in effect introducing some sort of 



7 RECOMMENDATIONS FOR FURTHER RESEARCH 



74 



quantum annealing. The relationship in Equation 12 that gives a sufficient 
evolution time so that the system is evolved adiabatically shows that the 
time is dependent on the minimum gap between the two lowest eigenvalues. 
By doubling the magnitude of the initial Hamiltonian, the gap at the start 
of the simulation is also doubled, and in many situations the initial gap is 
in fact the minimum gap over the whole simulation. With this increase of 
initial eigenvalue gap, the simulations in Table [9] have been re-run to yield 
the new results in Table HU 

In almost all cases the new probabilities exceed those of previous simulations. 
Indicating that greater gaps in eigenvalues do indeed reduce the possibility 
of Landau-Zener tunnelling to excited states. The system with no applied 
local Hamiltonian on the qubits has a minimum gap of 4 due to the initial 
Hamiltonian but increases towards the end of the simulation, so the prob- 
ability of measuring a correct state does in fact increase dramatically. The 
system where the local Hamiltonian has a magnitude of 1 on all qubits, the 
gap is 4 at both start and end of the simulation, as the gap does not in- 
crease more towards the end of the evolution, the probability of measuring 
a correct final state does increase but not as much as in the scenario of no 
local Hamiltonian. The simulation where the local Hamiltonian applied to 
all qubits has a magnitude of 2, still has a minimum gap of 2 due to the final 
Hamiltonian, although as the simulation started with a larger gap, the prob- 
ability of measuring a correct final state still increases. This small exercise 





Hamiltonian Magnitude 


Time 





1 


2 


0.1 


0.0080 (0.0082) 


0.0084 (0.0082) 


0.0084 (0.0082) 


1 


0.1337 (0.0608) 


0.1063 (0.0522) 


0.0543 (0.0342) 


2 


0.4245 (0.233) 


0.3228 (0.2024) 


0.1002 (0.0958) 


10 


0.9458 (0.8374) 


0.9334 (0.8496) 


0.9455 (0.8704) 


20 


0.9848 (0.9678) 


0.9822 (0.9698) 


0.9842 (0.9656) 



Table 11: Combined Final Probability of Correct Two States for Varying 
Time and Hamiltonian Magnitude for a Cube of Qubits with the Initial 
Hamiltonian Doubled in Magnitude where Brackets Indicate Previous Values 
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has shown that the eigenvalue gap between correct and incorrect states does 
affect the system undergoing adiabatic evolution, although it does seem that 
the dependence of total run time is not solely on the minimum gap, but the 
evolution of the gap. 

Therefore further research needs to be conducted on the possibilities of in- 
creasing Hamiltonian Magnitudes to yield greater eigenvalue gaps. 



7.2 Variable Rate Evolution 



Similarly to Section 7.1, due to probabilities increasing (although minimum 
gap being static in various scenarios) an investigation into the use of variable 
rate evolution may yield greater speed-up of computation, yet the question 
of at which point the speed-up should be applied (and for how long) during 
the evolution remains to be answered. 



7.3 Random Local Hamiltonian Magnitudes 



For the random 8 qubit simulation, there is a clear range of local Hamilto- 
nians that give extremely high probabilities for the correct states (viewable 



in Figure 25b). To induce eigenvalue gaps between correct states, using a 



randomly distributed set of local Hamiltonian within some finite range of 
magnitude on every qubit should prove effective. 



The random magnitudes for the various qubits were: 
hi = 1.78, h 2 = 1.8, h 3 = 1.85, h A = 1.73, h 5 = 1.77, h 6 



1.82, h 7 



1.9, h s 



1.69 



These were completely randomly selected, although care was taken to make 
sure the values did not differ much from the centre of Hamiltonian magnitude 
symmetry found in Section |5.5.2[ The values in brackets in Table [12] are 



those of previous bests taken from Table 10, therefore for T= 10, the values 
in brackets are previously from the local Hamiltonian magnitude of 1.75 and 
for T= 25 the previous magnitude being 2.25. It is clear that for T= 10, 
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State 


T= 10 


T=25 


111001101) 
10110011) 


0.1281 (0.2132) 
0.1701 (0.2132) 


0.1077 (0.2403) 
0.2034 (0.2403) 


101101101) 
01110011) 
10110110) 
11001110) 


0.1301 (0.1254) 
0.1375 (0.1254) 
0.1837 (0.1254) 
0.1881 (0.1254) 


0.0997 (0.1263) 
0.1176 (0.1263) 
0.2302 (0.1263) 
0.2212 (0.1263) 


Total 


0.9376 (0.9280) 


0.9798 (09858) 



Table 12: Table of Probabilities Obtained for All States Representing Max 
Independent Sets After Application of Randomly Distributed Local Hamil- 
tonian Magnitudes 

the probability of measuring a correct final state has in fact increased to 
a Figure better than that of the previous maximum, although for T= 25, 
the value has slightly decreased, even though it is itself an impeccably high 
probability. 

This study could be taken further to variations of the coupling strengths 
between qubits. Furthermore, this effect could be used to model errors in 
control parameters, as with all devices there is a finite internal error. 



7.4 More Highly Connected Graphs 

Looking at systems that are more highly connected than the standard cubic 
type is an option. This could lead to the development of a general algorithm 
across all types of graph for the MIS problem. By simply repeating the 



analysis in Section |4.3| for highly connected graphs, a simple mathematical 
proof may be deduced. This has been performed for all combinations in 
a system where all qubits have four interactions with nearest neighbours. 
The inequality developed to find the maximum Hamiltonian magnitude is: 
— 3H + 4 J < —5H — 4J. This gives an upper bound of H < 4. Study of 
the previous upper bound inequality for the cubic system along with this 
inequality gives a general upper bound inequality for all systems as: — (n — 
1)H + nJ < — (n + 1)H — nJ, where n is the number of nearest neighbour 
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interactions on each qubit. Therefore a general upper bound is given by H < 
n. Although a lower bound inequality must be developed and simulations of 
highly connected systems are appropriate. 

This research may also extend to the implementation of highly connected 
graphs onto a two dimensional planar chip. Optimisation algorithms do 
exist for the drawing of cubic planar graphs on regular lattices of qubits (see 
reference |2H]), although very little research has been conducted on problems 
where qubits are required to have a greater number of couplings. With chips 
that have been manufactured to date within the industry, the maximum 
number of couplings achieved is four, with edge and corner qubits having 
less. If this trend is set to continue, research into the some form of algorithm 
to implement problems on currently manufactured chips where qubits require 
more than four couplings is necessary. This will inevitably involve the use 
of strong ferromagnetic Ising interactions between neighbouring qubits that 
will effectively cause two qubits to act as one, with their combined remaining 
couplings accessible for needed problem interactions. 

7.5 Other NP-Hard Problems 

There are many variations of NP-Hard problems, and new problems can 
always be devised (whether there is any use for these newly devised prob- 
lems is another matter). However extensive lists are easily accessibk^] and 
algorithms for these other problems may be devised. 

A good example problem is an adaptation of the Exact Cover problem. Given 
a number of qubits, each of which can take a value or 1, and a set of clauses 
that each involves three qubits, whereby one qubit in every clause is of oppo- 
site orientation to the other two. Satisfying every clause will be an example 
of the adapted Exact Cover. Without researching the problem further and 
by carefully selecting the clauses so that no two clauses overlap by more than 
one qubit. The problem appears to be analogous to finding the ground state 



5 For example: 



http://www.csc.kth.se/~viggo/wwwcompendium/node276.html 
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Figure 26: 8 Qubits in a Randomly Selected Formation to Demonstrate an 
Adaptation of the Exact Cover Problem 

of a frustrated system, by connecting sets of three qubit in each clause into 
triangles and applying antiferromagnetic Ising interactions on every coupling, 
finding the ground state will give an adapted Exact Cover. 

As 8 qubits is the maximum that can be simulated, the following clauses 
of three qubits were randomly selected: 123, 346, 245, 158 and 678. These 
clauses allow for the system to be drawn in planar form, although it is im- 



portant to note that this system is not cubic, as seen in Figure 26 



After coding the system and running the simulation with an Q of 0.1, there 
were many states that ended the simulation with probabilities close to that 
of the maximum. After checking a selection of the highest probability state it 
was clear that they all satisfied all encoded clauses. A selection of those are 
as follows: |11010110), |11010101), |11010100), |11010011), |11010011) , etc. 
By including further clauses less states will end the simulation with high 
probability. 

This is quite a crude example, and only shown to work on a very specific 
system. By researching this further on a variety of systems, using a vari- 
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ety of clauses, only then can a general algorithm for this problem be fully 
developed. 
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8 Appendix 

8.1 Programming Code 
8.1.1 Instructions to User 

This code requires a Matlab installation of version 7.4 or greater. The file 
'SystemSolver.m' should be saved to a desired directory and this directory 
must be selected as the 'current directory' in the Matlab window. 

Initially the user must input data into the M-file depending on the type 
of system. The areas for this are designated in the code with examples of 
typical inputs. Running the M-file can be performed in one of two ways, if 
the M-file is open, the 'run' button can be pressed in this window. Otherwise 
typing 'run SystemSolver' in the main Matlab window will perform the same 
operation. 

After running the the user will be posed a series of questions. Firstly, a value 
of Q must be input by giving a typical number and pressing the 'enter' key. 
If no number is given, the code will use a default value of 0.1. Secondly the 
number of qubits in the system is input, this is obviously also dependent on 
the user and input in the same fashion. Finally, the simulation will ask the 
user to indicate the magnitude of the local Hamiltonian acting on the qubits, 
in most cases this will be the same for all qubits and can be input in this 
manner. Otherwise the user can indicate the individual local Hamiltonians 
before running the simulation in the M-file itself, the area for this is also 
clearly labelled with examples. 

Once theses inputs have been successfully completed the solver will continue 
to solve the system and output three to four plots and all final probabilities 
with corresponding states in the main Matlab window. 

There are two known bugs that may randomly crop up. Firstly, the pro- 
gram writes the system of ODEs to a separate M-File (named Hamilto- 
nian.ODE.m) in the current directory. If this file fails to be completed before 
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the program goes on the solve the system of ODEs stored in that file, the 
simulation will be interrupted. Running the simulation a second time will 
bypass this issue. Secondly, to plot all of the probabilities as a function of 
time, the plotting function is contained in a loop. The colour of each state's 
probability is determined by three random number generators giving indi- 
vidual values for Red, Green and Blue. An error may occasionally crop up 
at this point if two or more lines have identical colours, although this rests 
solely on all three random number generator generating identical numbers 
for all three inputs, all of which have three decimal places. The likelihood of 
this occurring is minimal although in some instances it could happen. 
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